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Abstract 

We illustrate a technique for specifying piecewise constant controls for classes of switched electrical 
networks, typically used in converting power in a dc-dc converter. This procedure makes use of decompo- 
sitions of SU{2) to obtain controls that are piecewise constant and can be constrained to be bang-bang 
with values or 1. Complete results are presented for a third order network first. An example, which 
shows that the basic strategy is viable for fourth order circuits, is also given. The former evolves on SO{3), 
while the latter evolves on 50(4). Since the former group is intimately related to 5(7(2) while the latter 
is related to SU{2) x SU{2), the methodology of this paper uses factorizations of SU{2). The systems 
in this paper are single input systems with drift. In this paper, no approximations or other artifices are 
used to remove the drift. Instead, the drift is important in the determination of the controls. Periodicity 
arguments are rarely used. 

Keywords: bang-bang controls, piecewise constant controls. Lie group, bilinear system, switched electri- 
cal network. 



1 Introduction 



In this paper the problem of explicit control of a class of switched electrical, lossless networks is considered. 
Specifically, it is shown how to determine explicitly piecewise controls, which can be constrained to take 
only the values 1 or 0, to achieve state transfers. Complete results are obtained for a third order lossless 
network, which has been studied before in jTl], ||^, |Q. The thesis, provides the model and assesses the 
controllability of the network. The paper, uses averaging to provide periodic controls for approximate 
state preparation. The same reference also emphasizes the desirability of finding hang-hang controls (with 
values 1 or Q), since this mode of control is closer to physical reality. In this paper a constructive protocol 
for precisely such a bang-bang control is provided. A fourth order network is also studied and preliminary 
results on certain explicit state transfers via bang-bang controls are provided. 

The energy conservation of the networks implies that they evolve on 5*0(3) (respectively 50(4)). For 
the third order network, the problem of bang-bang controls is susceptible to Euler factorizations (though 
non-Euler factorizations are also pertinent). However for constructiveness, explicit formulae, providing the 
Euler angles as expressions in the entries of the target state in 50(3), have to be supplied. To the best of our 
knowledge such explicit formulae are missing in the literature, especially when the two generators of so(3) (the 
Lie algebra of 50(3)), desired in the factorization, are the ones relevant to the model. It is worth emphasizing 
that the desired state in SO(3) does not already come specified with its Euler angles. Rather, it is described 
by the nine real entries which constitute this matrix. Similar issues (with the technicalities compounded) 
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present themselves for the fourth order network. It is primarily for this reason that the methodology of this 
paper uses a passage to an associated system evolving on SU{2) (respectively SU{2) x SU (2)). For the system 
associated to the third order network it turns out that Euler angles for SU{2), when the two generators are 
iax and iay, are needed. These are easier than the corresponding SO{3) angles to calculate because of two 
reasons: i) first, SU{2) matrices are 2x2 (the special unitarity mitigates the fact that the entries are complex) 
and thus, the matrix manipulations (which are inevitable if explicit formulae are required) are easier; and ii) 
SU{2) matrices admit the following representation (the Cayley-Klein representation): 

/ e^cosa e^'^sina \ . 
^ = ^("'^'^)=U^^^-"^«ina e-<cosa j- 

One such representation is nothing more than the entries written in polar coordinates. The advantage of 
(|l|) is that the condition SS* = S*S = I, det (S) — 1, is already incorporated. In contrast to SO{3), side 
conditions need not be stipulated. The attendant formulae, for even the SU{2) Euler angles, are messier if 
Cartesian coordinates were to be used (in our opinion, this is one of the reasons why explicit formulae for {x, y) 
Euler angles for SO{Z) are not available - there is no polar representation for real numbers). Furthermore, 
representing the columns of an £"0(3) matrix in spherical coordinates is equally unilluminating. In addition, 
for finding non-Euler factorizations, SU{2) is easier to work with. 

The differences between the third order network and the fourth order network examples are primarily 
twofold: i) for the fourth order network, factorizations of SU{2), different from ax,(Jy Euler angles, are 
needed. Indeed, the required factorizations are not of the Euler type. Such factorizations are easier to find 
when working with SU{2). ii) More importantly, the fourth order network problem amounts to the difficult 
question of constructive control of two systems with a single control. Due to the latter problem our results for 
the fourth order circuit are, pending further investigation, applicable under certain conditions on the circuit. 
Specifically, the transfers are achieved if any one of a set of relations between the constants of the circuits are 
satisfied. In part, these relations are a by-product of the specific choice of factorizations used. It should be 
possible to achieve these relations in practice, since they are only restrictions on the capacitors and inductors 
in the circuit. Work is ongoing to enlarge the class of state transfers and also to eliminate the restrictions on 
the constants. These preliminary results are, to the best of our knowledge, the first instances of constructive 
controllability for single input systems with drift evolving on SO (4). It is our opinion that, regardless of 
the specific model or the control technique, the most elegant manner to control a system evolving on 5*0(4) 
would indeed be to pass to an associated system on SU{2) x SU{2). Readers who are skeptical should first 
attempt to calculate the exponential of an so(4) matrix without any usage of SU (2) whatsoever. At a bare 
minimum manipulation of 4 x 4 matrices is required, whereas passage to SU (2) x SU{2) obviates all matrix 
manipulations. More importantly, finding e^, A £ so(4), via eigenvalues etc., occludes the structure of A in 
e^. This structure is relevant to the problem. 

Thus, the close relation between SU{2), 5*0(3) and 50(4) is used for the network systems. The group, 
SU{2), plays a prominent role in the control of many quantum systems (atoms and molecules. Cooper pairs, 
spin systems, photons and cxcitons). This explains the title of the paper. The rich algebraic structure of the 
Pauli matrices makes the deduction of the formulae easier than on the orthogonal groups. However, once a 
formula has been found on SU{2) - whether it be for an exponential or bang-bang controls etc., - it can be 
transferred easily to the orthogonal group. This is the rationale behind our method. 

Systems such as dc-dc switchmode power converters, in which switched electrical networks have a significant 
part, can be implemented in communication and data handling systems, portable battery-operated equipment 
and other applications. Thus, the results of this paper have useful consequences for these applications. Other 
strategies for controlling switched electrical networks use state-space averaging. Leonard and Krishnaprasad 
transform these systems into drift free systems and then apply averaging theory on Lie groups to specify 
small amplitude, periodic, open-loop controls for approximate state transfers. The approach of Sira-Ramirez 
, based on variable structure systems theory and sliding regimes, provides feedback controls for switched 
electrical networks. In contrast, the method in this paper obtains piecewise constant controls which further 
can be taken to be or 1, corresponding to the position of the switch. From the results of Jurdjevic and 
Sussmann it is known that bang-bang controls with values of and 1 can be used to prepare any target. 
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Thus, the paper provides constructive illustrations of the work in Q. It is emphasized that the approach 
taken in this paper does not resort to techniques for driftless systems by either i.) removing the drift via 
approximations or other methods which work only in fortuitous situations or ii.) by making use of periodicity. 
Arguments relying on periodicity are invalid in general ||^ and can lead to expensive controls even when valid. 
In this paper the only time periodicity is used is to rewrite free evolution terms with negative drift coefficients 
as free evolution terms with positive drift coefficients. 

The balance of this paper is organized as follows. In section 2, the relations between the unitary and 
orthogonal groups are reviewed. In the next section, the precise model for the 5*0(3) network is presented. 
Controls for this system are obtained in section 4. This section also contains the relevant formulae for the 
desired Euler angles. These are used to provide, first, piecewise constant controls and then bang-bang controls. 
The fifth section provides an illustration of the techniques for the fourth order network. The final section 
offers some conclusions. 



2 SU{2), 50(3) and 50(4) 

The Lie algebras su(2) = {V C^^^ | V* = -V} and so(3) = {VF G M^xs | y^T ^ „p^| ^^^^ isomorphic via 
the following explicit isomorphism, Q: 




(2) 



with 



1 
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-i 
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1 

-1 



Similarly, there is a group homomorphism (\) : SU{2) 50(3), [Q, obtained by considering the linear (vector 
space) map, Rjj : su{2) — > su{2), which for a fixed U E SU{2) is given by Rij{A) — UAU*. Identifying su{2) 
with i?'^, it can be shown that Ru G 50(3). The group homomorphism, (p, just associates U to Rjj- Finally, 
it can be shown, via a direct calculation using the Rodrigues' formula, that (/)(e^) — e'^'^^\ 

The groups 50(4) and SU{2) are related as follows. First, identifying the quaternions with via 
ei = 1, 62 = «, 63 = J, 64 = fc, leads to the following association, /, between a pair of unit quaternions, p, q 
and a linear map from to R^, [0 : 



/(p, q) — the linear map, 



pxq 



It can be shown that I{p,q) is an element of 50(4). Further, it is well known, Q, that the group of unit 
quaternions is explicitly isomorphic to SU{2). This then leads to a group homomorphism, : SU (2) x SU{2) 
50(4). It can be shown, via a direct calculation, that there is an associated Lie algebra isomorphism, 
ijj : su{2) X su{2) —> so(4), which satisfies (j){e^\e^'^) = e'^^^i^-^^)^ fo^. ^ny Ki x K2 in su{2) x su{2). -0 is 
given by 



4, {Ki X K2) 



( ' 
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bi 








(3) 



where 



§(ai+6i) 
(a2 + 62) + » (03 + 63)] 



H(«2 + &2) + j(a3 + &3)] 

-Uai+bi) 



(4) 
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and 



_ f |(6i-ai) ^ [(&2 - 02) + i (63 - aa)] \ 

^ i[-(fo2-a2) + z(63-a3)] -§(6i-ai) j ' ^ ^ 

Thus, given a system F = Ay + BVu{t), V G 5*0(3), one can associate a system, U = AU + BUu{t), 
U € SU{2), where A = ip~^i^) and S = to it. Now preparing a target, S in 5*0(3) with piecewise 

constant controls amounts to factoring S as ii'j^^_^Qi°-kA+bkB) ^ with > 0, if bk ^ 0. The condition, is 
either or 6^ = a^, is equivalent to preparing 5* with controls only taking values 1 or 0. As mentioned in 
the introduction, obtaining such factorizations explicitly is easier for SU{2). Hence, we work with the second 
system and factorize any matrix T in SU{2), such that (/)(T) = 5, as T = Hfc^i e("'=^+'"'^) with cither > 0, 
if bk 7^ etc. Recapitulating the preparation of a target 5 in 50(3) by associating it to a target T in SU{2), 



5 = 0(T)==</)^[|e'"=^+''^^j , 



(6) 



because is a homomorphism this gives 

Q 

5 = J] (e"-^+^'=^) , (7) 



fc=i 

and since Me^) = e'^(^) 



■Qga.A+fcfcS^ (8) 



k=l 

Therefore the same controls that prepare T also prepare S. The corresponding control values are, of course, 

bk. 
ak' 

Likewise, given a system V — AV + BVu{t), V G 50(4), two systems controlled by a single control 
u(t), are associated to it via, Ui = AiUi + BiUiu{t), Ui € 5*1/(2) and U2 = A2U2 + B2U2u{t), U2 G 51/(2). 
Here, (^1,^2) = -0"^^) and (-81,^2) = 0"^(i?). Given a target, 5* G 50(4), we prepare any (Ti,T2) such 
that (j)(Ti,T2) = 5. Usage of piecewise constant controls means that both the Ti have to be factorized as 
Ti = ll'j^^^eS"-''^^^^''^^\ with the same Q and same Ok and b^ for all fc = 1, . . . ,Q. The usual stipulations, 
afc > if 6fc 7^ etc., apply here too. 

Remark 1: It is well known, that the kernel of is {+I2, —I2} and that the kernel of is {(/2, 12) , {^h, ^^2)} ■ 
In this paper we do not make systematic use of this extra degree of freedom. 



3 The Third Order Network 



In this section, a switched electrical network with three circuit elements and no external constant power 
sources is considered. The switched electrical network examined here is identical to the one used by Leonard, 
Krishnaprasad and Wood [0,jn|. This network consists of two capacitors Ci and O2 with corresponding 
voltages Vi and V2 (see Figure |l ) . These capacitors are connected by a switch and an inductor L3 with current 
I3. The position of the switch connected to a control takes only the values of or 1. The control objective is 
to transfer energy from Ci to O2 by means of the inductor. Such systems, in the absence of external loads, 
can be modeled by defining the network state vector x = {xi,X2,X3)^ as xi = \fC\V\, X2 = VC2V2, and 
X3 = ^/LsIs- Let (jJi — Xj \JC\L3 and 102 = Xj^fC^L?,. Then the system is 



dt 






-wi(l 



u) 




[a + Bu) X. 



(9) 
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Figure 1: Third Order Switched Network 



If the control takes a constant vahie u for a time t, the state of the system can be written as 

x(t) = e(^+^")*x(0). 
The system on SU (2) associated to the system on 5*0 (3) , from section 2, is 



dt 



1 1 

--l^lO-y + - (iL!l(Jy + LlI2<7'x) U ] U. 



(10) 



(11) 



4 Quantum Control Techniques 

Preparing the final state, x^, is equivalent to the preparation of one of an infinite family of 50(3) matrices 
such that X/ = Sx (0) . This, in turn, defines a family of targets in SU{2) that can be associated with each 
such S. Each target S G S0{3) corresponds to two targets in SU{2). Preparing either of these two targets 
amounts to preparing S. 

A target T € SU{2), can be written as 



T — exp 



2 ("^^ 



bay + caz) 



(12) 



where a,b,c € R. Let A = ^Va'^ + 6^ + c^, s = — | {aax + bay + ca^) and p = ( a b c ) . An expression 
for the Lie group homomorphism, (j), described in section 2 is obtained from Rodrigues' formula |5|, 

1 — cos 2A rp 



r ; / \i T sin 2 A , , 

exp [tp (s)\ = 1 cos 2A H — — (s) 

2A 



4A2 



-PP 



(13) 



by using equation ^ and the fact that (p (T) = exp [ip (s)] . Now the problem of finding controls to drive the 
switched network on 5*0(3) from x (0) to x/ is converted to finding controls that prepare T £ SU{2). This is 
accomplished by using an appropriate decomposition of T. 



4.1 Decompositions of a target in SU(2) 

The general problem of preparing targets in SU{2) was considered in j^. In that work the theory requires 
that the drift and control matrices A and B be orthonormal. Orthonormality of A and B can be achieved by 
preliminary controls. The use of preliminary controls precludes the construction of bang-bang controls. Since 
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preliminary orthonormalization of the matrices A and B is not used here, the results do not follow directly 
from 1^. However, the general framework of that paper is helpful in this work. 

Consider the general problem of preparing a target for the system ( pi] ) in SU(2). The decompositions of 
elements of SU{2) considered in this paper are based on the fact that A and B in ( pi] ) are linear combinations 
of ia^ and iay. By writing the entries of T in the Cayley-Klein representation (|l|), various decompositions can 
be obtained. We describe three of these factorizations, one of which is used for general piecewise constant 
controls and the other two for bang-bang controls. First as shown in matrices in SU (2) may be decomposed 
into the following form: 

T{aX,^^)=e'P'''V{J)e'<•i-P^''', (14) 

for any p G R and 

/ ij \ 

V{-f)^exp[ ._ ' = exp [(-Im7) icTj^ + (Re7) icr:r] (15) 



17 

for 7 € C. In equation (|lj), 7 = a exp i (C + /i — 2p — ^) and p can be chosen so that V (7) is a free evolution 
factor. For the switched electrical network considered in this paper, the first and third factors of equation 
( p^ ) have a useful decomposition. In it was proved that the exponential of the third Pauli matrix can be 
expressed as a product of two factors 

6^^-^=^(71)1^(72), (16) 

where L e R. Thus it follows from equation (|l|) that T = 11^=1^(7^) > where 1 < Q < 5. In (|l^), let 
7fc = f 6*^*=, k — 1,2, then it holds that L ^ 9i — 62 + tt. Since L can be taken as an element of [0, 27r), it 
follows that |L — 7r| < TT, L 7^ 0. Thus for L 7^ 0, 71 and 72 can be chosen so that \9i — ^2! < tt. This means 
that 71 and 72 may be selected to lie any open half-plane. In other words, one can ensure that Ofc > 0, Vfc, as 

V (7fc) = exp bk'^ia^ + {bk - a^) ^icr.y . (17) 

It follows that the half-plane of interest is Im7fc > — ^Re7fc. This decomposition of the target in SU{2) 
provides piecewise constant controls with no further restrictions. 

Remark 2: The utility of the decomposition ( p^ ) is the following. Together with equation (|lj) it provides a 
factorization of any T in SU (2) of the type given by equation (H), with generally lower values of \/a1 + b'^ 
for each k, than would a factorization provided by the Euler parametrization (eqtns (|l^)-(|2^) below). This 
can be seen by viewing each factor in both of these decompositions as a matrix, V (7) (cf. eqtn (p^)). 

The complex numbers, 7, in the factorization provided by equation ( p^ ) each have a radial coordinate 
at most ^, whereas the radial coordinates due to Euler factorizations could be as high as 2n. This causes 
the former factorization to yield, generally, lower values for (individual and cumulative) \/ + ^fc- Since 
Ofc represents the duration and bk the power (=durationx amplitude) of the fcth pulse, this suggests that 
equation (^6|) is preferable for the simultaneous minimization of these two competing constraints, as long as 
it is reasonable to use any piecewise constant control. 



Next consider bang-bang controls. For free evolution, Uk = 0, bk = and from equation (|17|), V {jk) = 
exp (— afe^icTj,) , or in other words, 6k = When Uk — 1, this means that Uk = bk, so that V {jk) = 
exp (bk^icTx) which corresponds to the phase 6k = 0- These facts lead to the consideration of the following 
decomposition of the exponential of the third Pauli matrix: 



e^La, ^^^i^ay^iLa^^-rfay^ (18) 



The first and third factors of equation (g^ are free evolution and the second factor is obtained with a control 
pulse of 1. In each factor of (|lj) the drift coefficient is a positive number. With this decomposition the target 
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T from equation (p^, with p chosen so that ^ (7) is a free evolution factor, is prepared by at most seven 
factors of which at most two are control pulses. Now we consider another decomposition of T from which T 
is prepared by at most three factors. 

A different bang-bang protocol is obtained by the following. It is shown in [|j that the target T can be 
expressed as 

T(a,C,M) = e*^'^==e*-^'"«e*-^'"% (19) 
where D, E and F are solutions to the relations 



cos(i?) = ± Y cos^ C cos^ a + sin /zsin a (20) 
sin C cos a 

MD-F) = ± (21) 
ysm (, cos^ a + cos"^ /i sm a 

, „ sin /i sin a 

smiD + F) = ± J" (22) 

V cos^ Q cos^ a + sm /i sm a 



Remark 3: The expressions for the Euler angles given by equations (|20|)-(|22) were obtained by an explicit 
matrix calculation. Indeed, it is known that an Euler angle factorization, with factors that are exponentials of 
iuz and iuy, exists with a maximum of three factors. The orthogonality of the pairs {iuz, i<7y) and {iux, iuy) 
suggests an obvious Lie algebra isomorphism of su (2) with itself. This suggests that it should be possible to 
find a factorization of the type in equation ( p^ ) . This matrix calculation is facilitated by an explicit expression 
for the exponential of an su (2) matrix (which, incidently, is easier to manipulate than the corresponding so (3) 
expression). Even though there is a natural geometric equivalence between the pairs of generators of su (2) , 
(ifJz, i<^y) and {iux, ioy) , the expressions for the Euler angles are not simple consequences of one another. It is 
routine to show that the {iuz, i<yy) Euler angles are linear in the Cayley-Klein coordinates, whereas, equations 
(p0|)-(p2|) demonstrate that the {iax,i(7y) Euler angles involve transcendental functions. 



In the decomposition of T in equation ( |l9| ) the second factor is free evolution and control pulses of 1 are 
used to obtain the first and third factors. This decomposition of T prepares the target with at most three 
factors with no more than two control pulses. In summary, we have: 

Algorithm 1 Piecewise constant controls. 

1. Given an initial state x(0) and a final state xj, choose the numbers a, 5, and c in equation (|l3|) so that 
the target on 5*0(3), S ~ exp [V' (s)] , satisfies — Sx (0) . 

2. Write the entries of S in polar coordinates as in equation (|l]). 

3. Choose p in equation (|lj) so that V (7) is a free evolution factor. 

4. For each of the first and third factors of equation ( [T^ use the decomposition in equation (p^. Select the 
phases of 71 = f e'^^ and 72 — ^e*^^ in equation ([iq ) so that L = 9i — 62 + n and the drift coefficients 
are positive numbers. 



Algorithm 2 Bang-bang controls I 

Steps one through three are the same as for piecewise constant controls. 

Step four: For each of the first and third factors of equation ( |l4[ ) use the decomposition in equation ( |l8| 

Algorithm 3 Bang-bang controls II 

Steps one and two are the same as for piecewise constant controls. 

Step three: Solve for D, E and F in the decomposition of T given by equation (|l9|). 
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0.1, C2 = 0.2 and L3 = 0.5, then A 



lay 



and 



4.2 Example 

For the switched network in Figure |l|, if Ci 

_B = I (2\/5ay + V^Ocjx) ■ Suppose the initial state vector x(0) = (1,0,0)"^ and the final state vector x/ = 
(0, —1, 0)-^. Intermediate points for the system to traverse may be specified such as x(il) = (l/-\/2, 0, — 1/\/2)^ 
for the first intermediate point and x(i2) = (0, — 1/\/2, — l/v^)-^ for the second intermediate point Q. Sup- 
pose that it is desired for the system to pass through the intermediate points x(tl) and x(t2). Then three 
targets Ti, T2 and T3 in SU{2) must be prepared so that x(<l) = (/) (Ti)x(0), x(t2) = (jj {T2) yi{tl) and 
X/ = (/.(T3)x(i2). 

The target Ti is determined by the the Lie group homomorphism 



( C0S2A+ ° ^^4^^^^ 

ab(l-cos 2A) . esin 2 A 
lA^ 2A 

h sin 2A 



ab(l— cos 2A) 



i2A 



ac{l — cos 2A) 
4A2 



COS 2A 

bc(l-cos2A) 



V Jy? 2A 4X2" 

and X (tl) = (Ti) x(0) which lead to the relation 



2A 

b^(l-cos2A) 
iA^ 
g sin 2A 
2A 



^ cos2Ai 



(1 — cos 2Ai ) 



4A? 



a\b\ (1 — cos 2Ai ) 
4A^ 

aici (1— cos 2 Ai ) 
4A? 



ci sin 2Ai 

2Ai 
h\ sin 2Ai 

2Ai 



ae(l-cos 2A) , b sin 2A 
4A^ 2A 

be(l-cos2A) _ asin2A 
4P 



cos 2A 



2A 

c^(l-cos2A) 
4A5 




(23) 



(24) 



Also, ^ (Ti) e 50(3) requires that det [0 (Ti)] = 1. Let ai = and ci = 0, then 2Ai = and cos2Ai = 
Setting bi = J gives 



■(Ti) = 



1 





1 









1 





1 





1 


%/2 


V2 



(25) 



Therefore, Ti — exp (— -lio-j^) satisfies equation ( p^ and det [0 (Ti)] = 1. 

Similarly, the choice of the target T2 must satisfy x(t2) = 0(r2)x(tl) and 4'{T2) E SO{3). By letting 
a2 = and 62 = 0, then 2A2 = \c2\ and sin2A2 = ±1. Choosing C2 = — ^ leads to 



'(T2 




(26) 



Hence, r2 = exp (f icz) is a suitable target in SU{2). 

A target T3 that meets the requirements xy = (/)(T3)x(t2) and (/)(T3) e ^©(B) is obtained by letting 
63 = and C3 = from which it follows that 2 A3 = |a3| and cos 03 H — ^ sin 03 = — 1. Set = ~j and 



this gives 



V2' 



1 











^ ^ 

-— — 

" \/2 \/2 



(27) 



So Ts = exp (^iax) is a suitable target. 



4.2.1 Piecewise constant controls 



Algorithm^ is applied to the preparation of Ti, T2 and T3. By equations ( [l5| ) and (|T^), the real and imaginary 
parts of 7fc, fc = 1, . . . , Q are linear combinations of at and bk- For the network example, each of the factors 
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XI 



Figure 2: States obtained via general piecewise constant controls 
of T in equation (P) can be represented by 

exp bk]J^iax + (bk - ak) Vbiay 



(28) 



for k — 1, . . . ,Q. From this relation and afc > 0, we find that for each k, Im'jk > —\/2Kejk- So our choice of 
7A; must be above the line Im7fc = —^/2Re^k- This means that 9k € (— 0.304l7r, 0.69597r) , which is in keeping 
with the fact that 6k = f corresponds to free evolution and Ok — Q when = 1. 



8\/5' 



For the preparation of Ti = exp { — ^iay) , note that it can be achieved by free evolution with tl 
Now consider preparing T2 — exp (ficr^) . It follows from equation ( p^ ) that 

r2 = y (71) 1/(72), (29) 



where ^k — k — 1,2, and the phases of 71 and 72 must satisfy 62 - 

O2 — Using equation (Pq), the following coefficients are obtained for T2: 



Choose 



k 




h 


1 


0.649 


0.917 


2 


0.269 


-0.380 



and 

o 



(30) 



These coefficients, when placed in equation (||), indeed drive the system in 50(3) from x(il) to x (t2) as 
shown in Figure |^. 

To prepare ~ exp {-^iax) , observe that it can be obtained with a control pulse of one. Thus, tf — t2 = 

In Figure ^ the lines indicate that the network system is taken from a point in MP to another along some 
undetermined path on the unit sphere. Because 62 is a negative number, U2 is negative and, thus, U2 is not a 
control pulse that represents the position of the switch. 



4.2.2 Bang-bang controls I 

Now algorithm ^ is applied to the preparation of Ti, T2 and T3. As stated previously, the target Ti can be 
prepared by free evolution and the target T3 can be obtained with a control pulse of 1. Therefore, it remains 
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Figure 3: States obtained with bang-bang controls from algorithm]^ 

to get bang-bang controls for T2 so that the controls represent the position of the switch. It follows from 
equation (nsh that 



T2 = exp (i^cTz^ = e"'T'^«e*^'^- 



(31) 



Each of the three factors of equation (^Tj) are of the form of ( p8| ) , so that we have the following coefficients 
forTa: 



k 




bk 


1 
2 
3 


7r\/5 

20 
TvVTf) 

20 

20 




ir\/T0 
20 





(32) 



These coefficients represent bang-bang controls that drive the system as shown in Figure 

As in Figure H, each line in Figure ^ represents the system being taken from one point in M.^ to another 
along an undetermined path on the unit sphere. Controls having values of and 1 to prepare T2 can also be 
obtained by utilizing algorithm ^. 



4.2.3 Bang-bang controls II 

The preparation of Ti, T2 and is now achieved by use of algorithm ^ to obtain controls of and 1. Again, 
the target Ti can be prepared by free evolution and the target T3 can be obtained with a control pulse of 
1. It remains to prepare T2 — exp (i-|(Tz) by solving for D, E and F in the decomposition of T given by 
equation([ig|). Equation (|l|) is used to find values for a, C and 11. Since 



To = 



exp [ij) 






exp 



(33) 
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Figure 4: States obtained with bang-bang controls from algorithm |^ 
J 



V2 



Figure 5: Fourth Order Switched Network 



a = 0, C = f and /i G [0, 2tt) . From equations (po|), (|2^) and (|2^ ) we can choose D = 
F — ^. These values of D, E and F correspond to the following coefficients of T2: 



3ZL 

4 ' 



E = - 



k 


ak 




1 

2 
3 


3\/T07r 
20 

20 

20 


3\/T07r 
20 



20 



^ and 



(34) 



The switched network system in S0{3) is driven from x (tl) to x (t2) with these coefficients as shown in 
Figure ^ Each line in Figure ^ represents the system being taken from one point in MP to another along an 
undetermined path on the unit sphere. 



5 A Fourth Order Network 

In this section a fourth order network (see Figure ^) taken from [|ll| is considered and it is shown how to effect 
certain state transfers via bang-bang controls. To the best of our knowledge these state transfers are the first 



11 



instances of explicit exact control of systems with drift on S^, the sphere in R**. The system's equations are 

/0-i^0 0\ /O 07\ 

v 00(50 

-/3^+ 0-5 00 

\0 /3 0/ \-7 00/ 



xu{t) — Ax + Bxu (t) 



(35) 



The coefficient matrices of the system belong to so(4) and thus the system evolves on the sphere, S^, in M^. 
The constants, (3, 7 and 5 are positive and are related to the inductances and capacitances of the elements 
of the circuits. Specifically, we have, 

1 „ 1 1 „ 1 



V L1C2 y/LsUi V L1C4 V L3C2 

Here Ci,C2 are the two capacitances and ^1,^3 are the two inductances in the circuit (See for specific 
details). The state vector x — {xi,X2,X3,X4) is defined as xi — VL7/i, X2 — VC2V2, X3 — ^/L^I^ and 
X4 = y/UiVi. To this system we can associate two systems whose unitary generators evolve on SU{2) by 
using equations (||) - (||) of section 2: 



Ui = 



OzUl 



A1U1+ BiUiu{t) 



a^Uiu{t), Ui in SU{2) 



(36) 



and 



U2 = 



i (^Y^) ^-t^2 + * (^^) '^=cU2uit), U2 in SUi2) 

A2U2+B2U2u{t). 



(37) 



Note that both systems are controlled by the same control, u{t). 

The strategy to control the network is as follows. Supposed it is desired to transfer the state from 
(1,0,0,0)"^ to a vector y in 5*^, then one first represents (1,0,0,0)"^ and y by the unit quaternions 1 and y. 
The next step is to find a pair of unit quaternions p, q such that pq^^ = V- To p and q there correspond 
matrices (denoted by p and q again) in SU{2) (see section 2). We then try to find a u{t) which will prepare, 
simultaneously, p for system ( p6[ ) and q for system (^7|) . In general there will be an infinite family p, q such 
that pq~^ — y. To avail of this, we represent p via equation (|l|) with a, /i floating. This will then determine 
q. The parameters (a, C, m) are then found by the requirement that the same u{t) prepare both p and q. The 
details of this strategy, of course, depend on the specific sequence of piecewise constant controls which are 
used to prepare a state for a given system on SU(2). Equivalently, they depend on the specific factorization 
of SU (2) being employed. 

We will now illustrate one such technique with y = (0,0, 1,0)"'". The SU{2) matrix corresponding to y is 



e 2 



Thus, pq 



2 ^ y _ 



It turns out that this strategy can be implemented if any one of the following 



conditions on the constants of the circuit holds: 



= 2fc + 1 = 



7 + (5 



/c = 1,2,, 



(38) 



(3 — V ' (5 — 7' 

The above conditions are, of course, artifices of the specific factorization of SU{2) that will be presently 
employed. Note that these conditions imply that C2 — C4 and that = ^3^- 

Representing p as S (^j, ^4k(k+i) ' ~ ifc^fc+i^ ) suffices. Indeed, one can now factorize p as: 



p = cxp 
• exp 



(2fc + l)7r . \ / (2A: + l)7r ^ 
^1 exp ( '—Ai 



2(fc+ l)(z/ + /3) 
(2fc + 1) (6fc + l)7r 
2kJVTW) 



\V2{i^ + /3y 



(2fc + l)7r ^ 
V2{S + j) ' 



(39) 
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Since q is determined by p it follows that it can be factorized as: 



( (6fc + l)7r ^ \ 

Using equation ( ^8| ) it follows that the coefficients of Ai match those of A2 and similarly the coefficients 
of Bi match those of B2 ■ Since these coefficients represent the duration and power of the pieces of the control 
it follows that the same control, u{t) prepares both p and q and thus achieves the desired state for the 
circuit. Furthermore, the controls are indeed bang-bang with values 1 or 0. This follows from equation ( |38[ ) 
which forces, = 1 (keeping in mind the relation of these constants to the inductances and capacitances). 

Several remarks are in order at this stage: 



i) fc = was omitted from (p8[) since it would be physically unreasonable; 



ii) Using similar ideas, it can be shown that under equation (|38|), an explicit pulse sequence can be found for 
state transfer from (1, 0, 0, 0)"^ to any of the following states: ±(0, 1, 0, 0)^, ±(0, 0, 1, 0)^, ±(0, 0, 0, 1)'^. 



fi-v 



= 2fc + 1 = 



7+i 

7 — (5 



, k = 1,2,, 



Furthermore, this holds also when equation ( |38| ) is modified to 
For this it is useful to note that the Cay ley-Klein representation (|l|), need not necessarily be the polar 
coordinates of the entries. Indeed, since e*'^ = e^"^^ = — 1, one can begin with polar coordinates and yet 
dispense with the restriction that a be in [0, -IJ. To illustrate this, consider transferring the state from 
(1, 0, 0, 0)-^ to (0, 0, 0, 1)-^. Suppose equation (BSh is replaced with 



2fc-|-l: 



7 + (5 



ii — V (5 — 7 

Now pq^^ = e*^ . So if p = S* (a, , then 



(2fc+l) 



A:=l,2,. 



q = S[a--,--^,^--C 



(41) 



(42) 



It suffices to choose p = S (^j, ''^''^^(fc'+i) ' '^^^ffipfnr^) ' ^^^^ P ^® factorized as 

{2k + l)37r 



exp 



• exp 



2{k + l){iy + P) 
{2k +1) {I0k + 2)n 
2k{,y + p) 



Ai ) exp 

Ai 



({2k + l)j^ ^ {2k + 1). 



(43) 



Similarly, q can be factorized as 



exp 



37r 



• exp 



2{k + l){/3-iy) 
{10k + 2)7r 



A2 exp 



-A2 + 



V2{(3-y) V2{j~S) 



2k {13 - v) 



A, 



(44) 



Thus the same control which prepares p does likewise for q. Once again, this control is a bang-bang 
control with values and 1 ; 
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iii) The state (0, 1,0,0)"^ can be prepared by free evolution without any conditions on the circuit. This is 
not evident from equation ( |3^) , without the calculation of an exponential. On the other hand, since 
(0,1,0,0)"^ is equivalent to e^^"". it follows with minimal fuss upon passage to SU{2), i.e., with no 
calculation whatsoever. Indeed, the matrices Ai and A2 from equations ( p6|) and ( pTj) are (different) 
multiples of jcr^. Hence it suffices to choose the targets p and q of systems (|36|) and (p7[) as p = e'^°'^ 
and q = e^i^~^^'^' , for some L determined by 

iy + /3 l3-u' 

This p and q, and hence (0, 1, 0, 0)"^, can obviously be prepared by free evolution. Note, this conclusion 
did not even require an su (2) exponential. 

iv) For final states other than those in ii), we believe the same idea is viable. The resultant equations 
for the Cayley-Klein parameters of p now are transcendental Q . Intuitively, it seems plausible that 
these equations can be solved since 5*0(4) acts transitively on the 3-sphere, S*^, with isotropy given 
by SO (3) . As SO (3) is nearly SU{2), it seems reasonable to expect success of the strategy of finding, 
parametrically, a suitable p in SU{2) to ensure that the same u{t) will prepare both p and g; 

v) The conditions imposed by equation ( p^ ) or (^) can be met in practice. Nevertheless, it would be useful 
to achieve state transfers without this restriction. For this, different factorizations of SU{2) need to be 
developed; 

vi) The factorizations in equations (^9[), (^), ( |4^ ) and ( |4^ ) are not Euler factorizations. Even though three 
factors appear in each of these expressions, it is clear that Ai is not orthogonal to a linear combination of 
Ai and B^, i — 1,2. Thus, these factorizations are not Euler factorizations. Without passage to SU{2), 
it seems formidable to find similar factorizations directly for SO (4) . 

6 Conclusions 

Lossless networks of the type studied in this paper are important in many applications. Therefore, a con- 
structive strategy for preparing desired states in such circuits is interesting. In this paper the novel technique 
of using factorizations of the special unitary group SU{2) was shown to be a viable mechanism for this issue. 
The key enabling factor is the rich algebraic structure of su{2). Finding similar formulae for so(3) and so(4) 
directly is harder. However, once a formula on SU(2) has been found - whether it be for exponentials, bang- 
bang controls etc., - it can be transferred with ease to the orthogonal group. This is the rationale behind 
our method. While the individual properties of the networks played an important role in the success of the 
methodology, the basic idea of using factorizations of SU (2) is a useful complement to other methods for deal- 
ing with systems evolving on the unitary and orthogonal matrices. Indeed, for the treatment of systems with 
drift, techniques based on decompositions of unitary groups appear to be more viable. It seems reasonable to 
expect that similar methods should work for other systems evolving on the orthogonal groups. For instance, 
so(6) is isomorphic to su(4). The latter Lie algebra plays an important role in quantum control. 
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